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Abstract 

We present approximation algorithms for maximum independent set of pseudo-disks 
in the plane, both in the weighted and unweighted cases. For the unweighted case, we 
prove that a local search algorithm yields a PTAS. For the weighted case, we suggest 
a novel rounding scheme based on an LP relaxation of the problem, which leads to a 
constant-factor approximation. 

Most previous algorithms for maximum independent set (in geometric settings) 
relied on packing arguments that are not applicable in this case. As such, the analysis 
of both algorithms requires some new combinatorial ideas, which we believe to be of 
independent interest. 

1 Introduction 

Let F = {fx, . . . , /„} be a set of n objects in the plane, with weights W\, w 2 , . . . ,w n > 0, 
respectively. In this paper, we are interested in the problem of finding an independent set of 
maximum weight. Here a set of objects is independent, if no pair of objects intersect. 

A natural approach to this problem is to build an intersection graph G = (V, E), where 
the objects form the vertices, and two objects are connected by an edge if they intersect, 
and weights are associated with the vertices. We want the maximum independent set in G. 
This is of course an NP-Complete problem, and it is known that no approximation factor 
is possible within | | 1— e for any e>0ifNP^ZPP |Has96] . In fact, even if the maximum 
degree of the graph is bounded by 3, no PTAS is possible in this case |BF99j . 

In geometric settings, better results are possible. If the objects are fat (e.g., disks and 
squares), PTASs are known. One approach [Cha03j IEJS05] relies on a hierarchical spatial 
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subdivision, such as a quadtree, combined with dynamic programming techniques |Aro98j ; it 
works even in the weighted case. Another approach [Cha03j relies on a recursive application 
of a nontrivial generalization of the planar separator theorem |LT79tlSW98j ; this approach is 
limited to the unweighted case. If the objects are not fat, only weaker results are known. For 
the problem of finding a maximum independent set of unweighted axis-parallel rectangles, 
an 0(log log ^-approximation algorithm was recently given by Chalermsook and Chuzhoy 
[CC09j . For line segments, a roughly 0( y^Opt )-approximation is known [AM06] ; recently, 
Fox and Pach |FPllj have improved the approximation factor to n £ for not only for line 
segments but curve segments that intersect a constant number of times. 

In this paper we are interested in the problem of finding a large independent set in a set 
of weighted or unweighted pseudo-disks. A set of objects is a collection of pseudo-disks, if 
the boundary of every pair of them intersects at most twice. This case is especially intriguing 
because previous techniques seem powerless: it is unclear how one can adapt the quadtree 
approach [Cha03l IEJS05] or the generalized separator approach |Cha03] for pseudo-disks. 

Even a constant-factor approximation in the unweighted case is not easy. Consider the 
most obvious greedy strategy for disks (or fat objects): select the object e F of the smallest 
radius, remove all objects that intersect fi from F, and repeat. This is already sufficient to 
yield a constant-factor approximation by a simple packing argument |EKNS00l lMBH + 95] . 
However, even this simplest algorithm breaks down for pseudo-disks — as pseudo-disks are 
defined "topologically" , how would one define the "smallest" pseudo-disk in a collection? 

Independent set via local search. Nevertheless, we are able to prove that a different 
strategy can yield a constant-factor approximation for unweighted pseudo-disks: local search. 
In the general settings, local search was used to get (roughly) a A/4 approximation to 
independent set, where A is the maximum degree in the graph, see [Hal98] for a survey. In 
the geometric settings, Agarwal and Mustafa |AM06l Lemma 4.2] had a proof that a local 
search algorithm gives a constant-factor approximation for the special case of pseudo-disks 
that are rectangles; their proof does not immediately work for arbitrary pseudo-disks. Our 
proof provides a generalization of their lemma. 

In fact, we are able to do more: we show that local search can actually yield a PTAS for 
unweighted pseudo-disks! This gives us by accident a new PTAS for the special case of disks 
and squares. Though the local-search algorithm is slower than the quadtree-based PTAS in 
these special cases [Cha03j . it has the advantage that it only requires the intersection graph 
as input, not its geometric realization; previously, an algorithm with this property was only 
known in further special cases, such as unit disks |NHK 05j. Our result uses the planar 
separator theorem, but in contrast to the separator-based method in |Cha03j . a standard 
version of the theorem suffices and is needed only in the analysis, not in the algorithm itself. 

Planar graphs are special cases of disk intersection graphs, and so our result applies. 
Of course, PTASs for planar graphs have been around for quite some time [LT791 IBak94j . 
but the fact that a simple local search algorithm already yields a PTAS for planar graphs is 
apparently not well known, if it was known at all. 

We can further show that the same local search algorithm gives a PTAS for independent 
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set for fat objects in any fixed dimension, reproving known results in |Cha03t IEJS05] . 
This strategy, unfortunately, works only in the unweighted case. 

Independent set via LP. It is easy to extract a large independent set from a sparse 
unweighted graph. For example, greedily, we can order the vertices from lowest to highest 
degree, and pick them one by one into the independent set, if none of its neighbors was already 
picked into the independent set. Let do be the average degree in G. Then a constant fraction 
of the vertices have degree 0(dc), and the selection of such a vertex can eliminate 0(dc) 
candidates. Thus, this yields an independent set of size Q(n/da). Alternatively, for better 
constants, we can order the vertices by a random permutation and do the same. Clearly, 
the probability of a vertex v to be included in the independent set is l/(d(v) + 1). An easy 
calculation leads to Turan's theorem, which states that any graph G has an independent set 
of size >n/(d G + l) jA"S00] . 

Now, our intersection graph G may not be sparse. We would like to "sparsify" it, so that 
the new intersection graph is sparse and the number of vertices is close to the size of the 
optimal solution. Interestingly, we show that this can be done by solving the LP relaxation 
of the independent set problem. The relaxation provides us with a fractional solution, where 
every object fi has value Xi G [0, 1] associated with it. Rounding this fractional solution into 
a feasible solution is not a trivial task, as no such scheme exists in the general case. To this 
end, we prove a technical lemma (see Lemma 14.11) that shows that the total sum of terms of 
the form XiXj, over pairs fifj that intersect is bounded by the boundary complexity of the 
union of £ objects of F, where £ is the size of the fractional solution. The proof contains a 
nice application of the standard Clarkson technique |CS89j . 

This lemma implies that on average, if we pick fi into our random set of objects, with 
probability Xi, then the resulting intersection graph would be sparse. This is by itself suf- 
ficient to get a constant-factor approximation for the unweighted case. For the weighted 
case, we follow a greedy approach: we examine the objects in a certain order (based on a 
quantity we call "resistance"), and choose an object with probability around Xi, on condition 
that it does not intersect any previously chosen object. We argue, for our particular order, 
that each object is indeed chosen with probability Q(xi). This leads to a constant-factor 
approximation for weighted pseudo-disks. 

Interestingly, the rounding scheme we use can be used in more general settings, when 
one tries to find an independent set that maximizes a submodular target function. See 
Section 14.61 for details. 

Linear union complexity. Our LP analysis works more generally for any class of objects 
with linear union complexity. We assume that the boundary of the union of any k of these 
objects has at most gk vertices, for some fixed g. For pseudo-disks, the boundary of the 
union is made out of at most 6n — 12 arcs, implying g = 6 in this case [KLPS86] . 

A family F of simply connected regions bounded by simple closed curves in general 
position in the plane is k-admissible (with k even) if for any pair fi, fj 6 F, we have: (i) 
fi \ fj and fj \ fi are connected, and (ii) their boundary intersect at most k times. Whitesides 
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and Zhao [WZ90] showed that the union of such n objects has at most 3kn — 6 arcs; that 
is, q = 3k. So, our LP analysis applies to this class of objects as well. For more results on 
union complexity, see the sermon by Agarwal et al. [APS 08J. 

Our local-search PTAS works more generally for unweighted admissible regions in the 
plane. For an arbitrary class of unweighted objects with linear union complexity in the 
plane, local search still yields a constant-factor approximation. 

Rectangles. LP relaxation has been used before, notably, in Chalermsook and Chuzhoy's 
recent breakthrough in the case of axis-parallel rectangles |CC09j , but their analysis is quite 
complicated. Although rectangles do not have linear union complexity in general, we ob- 
serve in Section that a variant of our approach can yield a readily accessible proof of a 
sub logarithmic O (log n/ log log n) approximation factor for rectangles, even in the weighted 
case, where previously only logarithmic approximation is known |AvKS98l IBDMROTj ICha04] 
(Chalermsook and Chuzhoy's result is better but currently is applicable only to unweighted 
rectangles). 

Discussion. Local search and LP relaxation are of course staples in the design of ap- 
proximation algorithms, but are not seen as often in computational geometry. Our main 
contribution lies in the fusion of these approaches with combinatorial geometric techniques. 

In a sense, one can view our results as complementary to the known results on approxi- 
mate geometric set cover by Bronnimann and Goodrich [BG95] and Clarkson and Varadara- 
jan |CV07j . They consider the problem of finding the minimum number of objects in F to 
cover a given point set. Their results imply a constant-factor approximation for families 
of objects with linear union complexity, for instance. One version of their approaches is 
indeed based on LP relaxation |Lon01l IERS05j . The "dual" hitting set problem is to find 
the minimum number of points to pierce a given set of objects. Bronnimann and Goodrich's 
result combined with a recent result of Pyrga and Ray |PR08j also implies a constant-factor 
approximation for pseudo-disks for this piercing problem. The piercing problem is actually 
the dual of the independent set problem (this time, we are referring to linear programming 
duality). We remark, however, that the rounding schemes for set cover and piercing are based 
on different combinatorial techniques, namely, e-nets, which are not sufficient to deal with 
independent set (one obvious difference is that independent set is a maximization problem). 

In Theorem 14.61 we point out a combinatorial consequence of our LP analysis: for any 
collection of unweighted pseudo-disks, the ratio of the size of the minimum piercing set to 
the size of maximum independent set is at most a constant. (It is easy to see that the ratio 
is always at least 1; for disks or fat objects, it is not difficult to obtain a constant upper 
bound by packing arguments.) This result is of independent interest; for example, getting 
tight bounds on the ratio for axis-parallel rectangles is a long-standing open problem. 

In an interesting independent development, Mustafa and Ray |MR09] have recently ap- 
plied the local search paradigm to obtain a PTAS for the geometric set cover problem for 
(unweighted) pseudo-disks and admissible regions. 



4 



2 Preliminaries 



In the following, we have a set F of n objects in the plane, such that the union complexity of 
any subset X C F is bounded by q \X\, where q is a constant. Here, the union complexity 
of X is the number of arcs on the boundary of the union of the objects of X. Let A(F) 
denote the arrangement of F, and V (F) denote the set of vertices of .A(F). 

In the following, we assume that deciding if two objects intersect takes constant time. 

3 Approximation by Local Search: Unweighted Case 

3.1 The algorithm 

In the unweighted case, we may assume that no object is fully contained in another. 

We say that a subset L of F is b-locally optimal if / is an independent set and one 
cannot obtain a larger independent set from / by deleting at most b objects and inserting at 
most 6+1 objects of F. 

Our algorithm for the unweighted case simply returns a Mocally optimal solution for a 
suitable constant b, by performing a local search. We start with L ^— 0. For for every subset 
X C F \ L of size at most b + 1, we verify that X by itself is independent, and furthermore, 
that the set Y C L of objects intersecting the objects of X, is of size at most |X| — 1. If 
so, we do L 4— (L \ Y) U X. Every such exchange increases the size of L by at least one, 
and as such it can happen at most n times. Naively, there are ( 6 "J subsets X to consider, 
and for each such subset X it takes 0(nb) time to compute Y. Therefore, the running time 
is bounded by 0(n 6+3 ). (The running time can be probably improved by being a bit more 
careful about the implementation.) 

3.2 Analysis 

We present two alternative ways to analyze this algorithm. The first approach uses only the 
fact that the union complexity is low. The second approach is more direct, and uses the 
property that the regions are admissible. 

3.2.1 Analysis using union complexity 

The following lemma by Afshani and Chan |AC06] , which was originally intended for different 
purposes, will turn out to be useful here (the proof exploits linearity of planar graphs and 
the Clarkson technique |CS89j ): 

Lemma 3.1 Suppose we have n disjoint simply connected regions in the plane and a col- 
lection of disjoint curves, where each curve intersects at most k regions. Call two curves 
equivalent if they intersect precisely the same subset of regions. Then the number of equiva- 
lent classes is at most c§nk 2 for some constant cq. 
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Let S be an optimal solution, and let L be a 6-locally optimal solution. We will upper- 
bound |S| in terms of |L|. 

Let S>6 denote the set of objects in S that intersect at least 6+1 objects of L. Let §<b 
be the set of remaining objects in S. 

If ft G S intersects fj G L, then the pair of objects contributes at least two vertices to 
the boundary of the union of S U L. Indeed, the objects of S (resp. L) are disjoint from 
each other since this is an independent set, and no object is contained inside another (by 
assumption). We remind the reader that for any subset ICF, the union complexity of the 
regions of X is < g \X\. As such, the union complexity of S>& U L is < g (|S>&| + |L|). Thus, 

2(6+l)|S> 6 | < e (|S> 6 | + |L|) 

Q n i 



»>6 



< 



On the other hand, by applying Lemma 13.11 with L as the regions and the boundaries of 
S<b as the curves, the objects in S<b form at most co6 2 |L| equivalent classes. Each equivalent 
class contains at most b objects: Otherwise we would be able to remove b objects from L 
and intersect b + 1 objects in this equivalence class to get an independent set larger than L. 
This would contradict the 6-local optimality of L. Thus, |S<&| < co6 3 |L|. 

Combining the two inequalities, we get 

|S| < |S< 6 | + |S >6 | < ( c b 3 + 



2(6+1) - g, 

For example, we can set b = \g/2] and the approximation factor is 0(g 3 ). 

Theorem 3.2 Given a set ofn unweighted objects in the plane with linear union complexity, 
for a sufficiently large constant b, any b-locally optimal independent set has size fi(opt), where 
opt is the size of the maximum independent set of the objects. 



3.2.2 Better analysis for admissible regions 

A set of regions F is admissible, if for any two regions /, /' G F, we have that f \ f and 
f'\f are both simply connected (i.e., connected and contains no holes). Note that we do 
not care how many times the boundaries of the two regions intersect, and furthermore, by 
definition, no region is contained inside another. 

Lemma 3.3 Let F be a set of admissible regions, and consider a independent set of regions 
ICF, and a region f G F \ /. Then, the core region f \ I = f \ U g is non-empty and 

simply connected. 

Proof: It is easy to verify that for the regions of I to split / into two connected components, 
they must intersect, which contradicts their disjointness. ■ 
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Lemma 3.4 Let X, Y C F be two independent sets of re- 
gions. Then the intersection graph G of X U Y is planar. 



Proof: Lemma 13.31 implies the planarity of this graph. 

Indeed, for a region / £ X, the core f' = f\ U g is 

non-empty and simply connected. Place a vertex Vf inside 
this region, and for every object g £ Y that intersects /, 
create a curve from Vf to a point p^ g on the boundary of 
g that lies inside /. Clearly, we can create these curves in 
such a way that they do not intersect each other. Pf9 ' , 

Similarly, for every region g £ Y, we place a vertex v g inside g, and connect it to all the 
points pf t9 placed on its boundary, by curves that are contained in g, and they are interior 
disjoint. Clearly, together, these vertices and curves form a planar drawing of G. ■ 

We need the following version of the planar separator theorem. Below, for a set of vertices 
U in a graph G, let T(U) denote the set of neighbors of U, and let T(U) = T(U) U U . 




Lemma 3.5 ( [Fre87| ) Th ere are constants c\, C2 and c 3 , such that for any planar graph 
G = (V, E) with n vertices, and a parameter r, one can find a set of X C V of size at most 
C\n/ y/r, and a partition of V \ X into n/r sets VI, ... , V n / r} satisfying: (i) \Vi\ < C2r, (ii) 
T(Vi) n Vj = 0, for i ^ j, and (in) \T(Vi) n X\ < c 3v ^. 

Let S be the optimal solution and L be a 6-locally optimal solution. Consider the bipartite 
intersection graph G of SUL. By Lemma l3~4"l we can apply Lemma l3~5l to G, for r = fr/( c 2+ c 3)- 
Note that |r(Vi)| < c 2 r + c 3v /r < b for each i. Let 

Si = \Vi n S| , 

ii = \Vi n L| , 

and bi = \T(Vi) fl X\ , for each i. 

Observe that ii + bi > Sj, for all %. Indeed, otherwise, we can throw away the vertices of 
L fl T(Vi) from L, and replace them by V{ fl S, resulting in a better solution. This would 
contradict the local optimality of L. Thus, 



|S| < ^Si + \X\ < 



\x\ 



< LI + c 3 vr • 



< 1 1_ I + (d + c 3 ; 



r 

ISI + 1 1_ I 



ex- 



it follows that |S| 
following. 



< 



'1 + 0(l/Vb))\L\. We can set b to the order of 1/e 2 , and we get the 
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Theorem 3.6 Given a set of n unweighted admissible regions in the plane, any b-locally 
optimal independent set has size > (1 — 0(l/\/b))opt, where opt is the size of the maximum 
independent set of the objects. In particular, one can compute an independent set of size 
> (1 — e)opt, in time n°( l l £ \ 



3.2.3 Analysis for fat objects in any fixed dimension 

We show that the same algorithm gives a PTAS for the case when the objects in F are fat. 
This result in fact holds in any fixed dimension d. For our purposes, we use the following 
definition of fatness: the objects in F are fat if for every axis-aligned hypercube B of side 
length r, we can find a constant number c of points such that every object that intersects B 
and has diameter at least r contains one of the chosen points. 

Smith and Wormald |SW98] proved a family of geometric separator theorems, one version 
of which will be useful for us and is stated below (see also |Cha03] ): 



Lemma 3.7 (| ISW9 8j) Given a collection of n fat objects in a fixed dimension d with con- 



stant maximum depth, there exists an axis-aligned hypercube B such that at most 2n/3 objects 
are inside B, at most 2n/3 objects are outside B, and at most 0(n l ~ l l d ) objects intersect 
the boundary of B. 

We need the following extension of Smith and Wormald's separator theorem to multiple 
clusters (whose proof is similar to the extension of the standard planar separator theorem in 
[Fre87] ): 

Lemma 3.8 There are constants c\, C2, C3 and C4, such that for any intersection graph 
G = (V, E) of n fat objects in a fixed dimension d with constant maximum depth, and a 
parameter r, one can find a set of X C V of size at most c\n/r l / d , and a partition ofV\X 
into n/r sets Vi, . . . , V n / r , satisfying: (i) \V,j\ < c 2 r , (ii) T(Vi) Pi Vj — 0, for i 7^ j , and (Hi) 
J2i |rrVi) n X\ < c 3 n/rV d ; and (iv) \T(Vi) H X\ < c 4 r. 

Proof: Assume that all objects are unmarked initially. We describe a recursive procedure 
for a given set S of objects. If \S\ < c 2 r, then S is a "leaf" subset and we stop the recursion. 
Otherwise, we apply Lemma I3T71 Let S' and 5"' be the subset of all objects inside and outside 
the separator hypercube B respectively. Let S be the subset of all objects intersecting the 
boundary of B. We mark the objects in S and recursively run the procedure for the subset 
S' U S and for the subset S" U S. 

Note that some objects may be marked more than once. Let X be the set of all objects 
that have been marked at least once. For each leaf subset Si, generate a subset Vi of all 
unmarked objects in Si. Property (i) is obvious. Properties (ii) and (iv) hold, because the 
unmarked objects in each leaf subset Si can only intersect objects within Si and cannot 
intersect unmarked objects in other S/s. 

The total number of marks satisfies the recurrence T{n) = if n < c 2 r, and 

T(n) < max T(m + 0(n 1_1 / d )) + T{n 2 + 0{n^ 1/d )) + 0{n^ lld ) 

ni ,n,2<2n/3 



S 



otherwise. The solution is T(n) = 0(n/r l l d ). Thus, we have \X\ = 0(n/r 1 ' d ). Furthermore, 
for each object / 6 X, the number of leaf subsets that / is in is equal to 1 plus the number 
of marks that / receives. Thus, (iii) follows. ■ 

Let S be the optimal solution and L be a 6-locally optimal solution. Consider the bipartite 
intersection graph G of S U L, which has maximum depth 2. We proceed as in the proof 
from Section 13.2. 2[ using Lemma 13.81 instead of Lemma 13.51 Note that (iii)-(iv) are weaker 
properties but are sufficient for the same proof to go through. The only differences are that 
square roots are replaced by cZth roots, and we now set r = b/(c 2 + c 4 ), so that |r(V^) < 
c 2 r + c 4 r < b. We conclude: 

Theorem 3.9 Given a set of n fat objects in a fixed dimension d, any b-locally optimal inde- 
pendent set has size > (1 — 0(l/6 1 / d ))opt , where opt is the size of the maximum independent 
set of the objects. In particular, one can compute an independent set of size > (1 — e)opt, 
in time n ol ^l e '. 



4 Approximation by LP Relaxation: Weighted Case 
4.1 The algorithm 

We are interested in computing a maximum-weight independent set of the objects in F = 
{/i, . . . , /„}, with weights Wi, . . . , w n , respectively. To this end, let us solve the following LP 
relaxation: 

n 

max WiXi (1) 
J>*<1 VpGV(F) 

fiBp 

< Xi < 1, 

where V (F) denotes the set of vertices of the arrangement .A(F). 

In the following, X{ will refer to the value assigned to this variable by the solution of the 
LP. Similarly, Opt = Yli w i x i wm denote the weight of the relaxed optimal solution, which 
is at least the weight opt of the optimal integral solution. 

We will assume, for the time being, that no two objects of F fully contain each other. 

For every object f i: let its resistance be the total sum of the values of the objects that 
intersect it. Formally, we have that 

v(fuF)= Yl x r 
fj£F\{fi} 

We pick the object in F with minimal resistance, and set it as the first element in the 
permutation II of the objects. We compute the permutation by performing this "extract- 
min" , with the variant that objects that are already in the permutation are ignored. Formally, 
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if we computed the first % objects in the permutation to be LL = (tti, . . . ,7Tj), then the (i+l)th 
object 7Tj +1 is the one realizing 

ri i+1 = min V (f, F\n i ). (2) 

The algorithm starts with an empty candidate set C and an empty independent set /, 
and scans the objects according to the permutation in reverse order. At the zth stage, the 
algorithm first decides whether to put the object 7r n _j in C, by flipping a coin that with 
probability x(7r n _j)/r comes up heads, where x(7r n _j) is the value assigned by the LP to the 
object 7r n „j and r is some parameter to be determined shortly. If 7r n _j is put into C then we 
further check whether 7r n _j intersects any of the objects already added to the independent 
set /. If it does not intersect any objects in J, then it adds ir n ^i to / and continues to the 
next iteration. 

In the end of the execution, the set / is returned as the desired independent set. 



4.2 Analysis 

Let F be a set of n objects in the plane, and let u(m) be the maximum union complexity of 
m < n objects of F. Furthermore, we assume that the function u(-) is a monotone increasing 
function which is well behaved; namely, u(n) /n is a non-decreasing function, and there exists 
a constant c, such that u(xr) < cu(r), for any r and 1 < x < 2. In the following, a vertex p 
of V(H) is denoted by (p, i,j), to indicate that it is the result of the intersection of the ith 
and jth object. 

The key to our analysis lies in the following inequality, which we prove by adapting the 
Clarkson technique |CS89j . 

Lemma 4.1 Let H be any subset of F. Then XiXj = 0(u(£(H))) ; where £(H) = 

( P ,ij)eV(H) 

Proof: Consider a random sample R' of H, where an object fa is being picked up with 
probability Xj/2. Clearly, we have that (p,i,j) e V(F) appears on the boundary of the 
union of the objects of R', if and only if fa and fj are being picked, and none of the objects 
that cover p are being chosen into R'. In particular, let U(R') denote the vertices on the 
boundary of the union of the objects of R'. We have 



Pr 



e U(R') 



X 3 



n 

fk3p, 



1 _ ^± 

2 



> 



y 



> 



by the inequality Y\ k iX ~ a k) > 1 — J2k a k f° r a k £ [0,1], since ^2f k3p Xk < 1 (as the LP 
solution is valid). On the other hand, the number of vertices on the union is |U(R')| < u(|R'|). 
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Thus, 



E 

(p,*,i)6V(F) 



00 iX j 



< Pr[(p,i,j')eU(R') 

(pAi)eV(F) 



E 



|U(R') 



< E 



u(|R'l) 



To bound last expression, observe that \i = E[|R'|] = £(H)/2. Furthermore, by Chernoff in- 
equality, Pr[|R'| > (t + l)(i] < 2~*. Thus, E[u(|R'|)] < E t =i 2"* +1 u(t/i) = Et=i 2- m t° (1) u( / u) 
0(u(/i)) = 0(u(£(H))), since u(-) is well behaved. ■ 

Lemma 4.2 For any i, the resistance of the ith object m (as defined by Eq. is rji = 
u(8(F))' 



O 



£(R 



Proof: Fix an z, and let K = F \ {7Ti, . . . , 7Tj_i}. By Lemma 14. 14 

X)x i i7(/ i) K) = 2 £ x iXj = 0(u(£(K))) 
/ 3 eK (p,*j)eV(K) 



E 



£(K) 



*7(/i,K) = 



u(g(K)) 
£(K) 



u(£(F)) 
8(F) 



by the monotonicity of u(n) /n, where £(K) = XI/ gK-^i- ^ follows that rnin/ j6 K vifji K) 

^ £(F) 

Lemma 4.3 For a sufficiently large constant c, setting r = c "^ , £/ie algorithm in Sec- 
tion 4-1 outputs in expectation an independent set of weight Q((n/u(n))Opt). 

Proof: Indeed, the j'th object in the permutation is added to C with probability xj/t. Let 
K be the set of all objects of F that were already considered and intersect Xj. Clearly, £(K) 
is exactly the resistance rjj of fj. Furthermore by picking c large enough, we have that 
£(K) = rjj < r/2 by Lemma [4.21 This implies that 



Pr 



fit* 



f)eC = Pr [k n C = fj g C = Pr[K n C = 0] = n( 1- ^ 



/,eK 



£ K 1 

1 ^ > -, 

r r - 2' 



by the inequality rifc(l ~ a k) > 1 _ a k fo r °fc e [0, !]■ Now, we have that 



% = Prft 6 /] = Pr [/, G / £ 6 C\ • Pr[/, GC]>^. 
As such, the expected value of the independent set output is 



as 8(F) < n. 
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4.3 Remarks 



Variant. In the conference version of this paper |CH09j . we proposed a different variant 
of the algorithm, where instead of ordering the objects by increasing resistance, we order 
the objects by decreasing weights. An advantage of the resistance-based algorithm is that 
it is oblivious to (i.e., does not look at) the input weights. This feature is shared, for 
example, by Varadarajan's recent algorithm for weighted geometric set cover via "quasi- 
random sampling" [VarlOj . Another advantage of the Resistance-based algorithm is its 
extendibility to other settings; see Section 14.61 

Derandomization. The variance of the expected weight of the returned independent set / 
could be high, but fortunately the algorithm can be derandomized by the standard method of 
conditional probabilities/expectations |MR95j . To this end, observe that the above analysis 
provide us with a constructive way to estimate the price of generated solution. It is now 
straightforward to decide for each region whether to include it or not inside the generated 
solution, using conditional probabilities. Indeed, for each object we compute the expected 
price of the solution if it is include in the solution, and if it is not included in the solution, 
and pick the one that has higher value. 

Coping with object containment. We have assumed that no object is fully contained 
in another, but this assumption can be removed by adding the constraint X]/ c/ x i — ^ ^ or 
each i to the LP. Then, for any subset H of F, we have 

^ x t Xj < £(H), 

fi,fj€H 
fiCfi 

and so Lemma 14.11 still holds. The rest of the analysis then holds verbatim. 

Time to solve the LP. This LP is a packing LP with 0(n 2 ) inequalities, and n variables. 
As such, it can be (1 + ^-approximated in 0(n 3 + e~ 2 n 2 log n) = 0(n 3 ) by a randomized 
algorithm that succeeds with high probability |KY07j . For our purposes, it is sufficient to 
set e to be a sufficient small constant, say e = 10~ 4 . 

We have thus proved: 

Theorem 4.4 Given a set of n weighted objects in the plane with union complexity 0(u(n)), 
one can compute an independent set of total weight f2((n/u(n))opt), where opt is the max- 
imum weight over all independent sets of the objects. The running time of the randomized 
algorithm is 0(n 3 ), and polynomial for the deterministic version. 

Corollary 4.5 Given a set of n weighted pseudo-disks in the plane, one can compute, in 
0(n 3 ) time, a constant factor approximation to the maximum weight independent set of 
pseudo-disks. 
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Theorem 14.41 can be applied to cases where the union complexity is low. Even in the case 
of fat objects, where PTASs are known [Cha03irEJS05j . the above approach is still interesting 
as it can be extended to more general settings, as noted in Section 14.61 

4.4 A combinatorial result: piercing number 

In the unweighted case, we obtain the following result as a byproduct: 

Theorem 4.6 Given a set ofn pseudo-disks in the plane, let opt be the size of the maximum 
independent set and let opt' be the size of the minimum set of points that pierce all the 
pseudo-disks. Then opt = fi(opt'). 

Proof: By the preceding analysis, we have opt = fi(Opt), i.e., the integrality gap of our LP 
is a constant. (Here, all the weights are equal to 1.) 
For piercing, the LP relaxation is 

n 

min ^2 yp 

pGV(S) 

- 1 Vz = l,...,n 

pe/i 

< y p < 1. 

Let Opt' be the value of this LP. Known analysis |Lon01irERS05] implies that the integrality 
gap of this LP is constant if there exist e-nets of linear size for a corresponding class of 
hypergraphs formed by objects in F and points in V (S). Pyrga and Ray [PR08t Theorem 12] 
obtained such an existence proof for this ("primal") hypergraph for pseudo-disks. Thus, 
opt' = O(Opt'). 

To conclude, observe that the two LPs are precisely the dual of each other, and so 
Opt = Opt'. ■ 

4.5 A discrete version of the independent set problem 

We now show that our algorithm can be extended to solve a variant of the independent set 
problem where we are given not only a set F of n weighted objects but also a set P of m 
points. The goal is to select a maximum-weight subset § C F such that each point p G P is 
contained in at most one object of S. (The original problem corresponds to the case where 
P is the entire plane.) Unlike in the original independent set problem, it is not clear if local 
search yields good approximation here, even in the unweighted case. 

We can use the same LP as in Section 14.11 to solve this problem, except that we now 
have a constraint for each p 6 P instead of each p 6 V(F). In the rest of the algorithm and 
analysis, we just reinterpret D fj ^ 0" to mean Pi fj D P ^ 0". 

Lemma 14.11 is now replaced by the following. 
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Lemma 4.7 Let H be any subset of F. T/ien sCjSCj = 0(u(£(H))). 

/i,/jeH 

Proof: Consider a random sample R of H, where an object f\ is being picked up with 
probability £;/2. Let VT>(R) be the cells in the vertical decomposition of the complement of 
the union 1X(R). For a cell A G VD(R), let xa = X^gh f 4 nmt(A)^0 x « be ^ ne total energy of the 
objects of H that intersects the interior of A. A minor modification of the analysis of Clarkson 
|CS89] implies that, for any constant c, it holds E ^AeVD(R) ( X ^-T = ^(e[u(£(H))] Y 

A point of p G P is active for R, if it is outside the union of objects of R. Let P' be the 
set of active points in P. We have that 



E 

/in/,nP'^0 



OC 2 OC j 



E 

AeVD(R) 



X, 



Furthermore, by arguing as in Lemma I4.1[ every point p G P has probability at least 
II/iGH peA^ ~~ x i/fy > 1/2 to be active. Thus, 



1 x ^ 

2 XiX 'j — ^ 




< E 


E 4 








AeVD(R) 



0(E[u(£(H))]) =0(u(£(H))) 



again, by arguing as in Lemma [4.11 



The above proof is inspired by a proof from |AACS 98j . There is an alternative argument 
based on shallow cuttings, but the known proof for the existence of such cuttings requires a 
more complicated sampling analysis |Mat92] . 

The rest of the analysis then goes through unchanged. We therefore obtain an 0(1)- 
approximation algorithm for the discrete independent set problem for unweighted or weighted 
pseudo-disks in the plane. 



4.6 Contention resolution and submodular functions 



The algorithm of Theorem 14.41 can be interpreted as a contention resolution scheme; see 
Chekuri et al. [CVZllj for details. The basic idea is that given a feasible fractional solution 
x G [0, l] n , a contention resolution scheme scales down every coordinate of x (by some 
constant b) such that given a random sample C of the objects according to x (i.e., the ith 
object fi is picked with probability bxj), the contention resolution scheme computes (in our 

case) an independent set I such that Pr fi & I fi G C > c, for some positive constant c. 

The proof of Lemma 14.31 implies exactly this property in our case. 

As such, we can apply the results of Chekuri et al. [CVZllj to our settings. In partic- 
ular, they show that one can obtain constant approximation to the optimal solution, when 
considering independence constraints and submodular target function. Intuitively, submod- 
ularity captures the diminishing-returns nature of many optimization problems. Formally, 
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a function g : 2 F -> ]R is submodular if #(X U F) + g(X n F) < g(X) + g(Y), for any 
X,F C F. 

As a concrete example, consider a situation where each object in F represents a coverage 
area by a single antenna. If a point is contained inside such an object, it is fully serviced. 
However, even if it is not contained in a object, it might get some reduced coverage from the 
closest object in the chosen set. In particular, let u(r) be some coverage function which is the 
amount of coverage a point gets if it is at distance r from the closest object in the current 
set /. We assume here v(-) is a monotone decreasing function. Because of interference 
between antennas we require that the regions these antennas represent do not intersect (i.e., 
the set of antennas chosen needs to be an independent set). 

Lemma 4.8 Let P be a set of points, and let F be a set of objects in the plane. Let u(-) be 
a monotone decreasing function. For a subset H C F, consider the target function 

a P (H) = J>(clp(H)) , 
peP 

where d p (H) is the distance of p to its nearest neighbor in H. Then the function ap(H) is 
submodular. 

Proof: The proof is not hard and is included for the sake of completeness. For a point p e P, 
it is sufficient to prove that the function z/(d p (H)) is submodular, as «p(H) is just the sum 
of these functions, and a sum of submodular functions is submodular. 

To prove the latter, it is sufficient to prove that for any sets X C F C F, and an object 

/eF\F, 

u(d p (X U {/})) - zz(dpPO) > ^(d P (F U {/})) - u(d p (Y)) . 

To this end, let x and y be the closest objects to p in X and F, respectively. Similarly, let 
£ x ,£ y ,£f be the distance of p to x,y and /, respectively. The above then becomes 

v(wm{£ x ,e f )) - 1/(4) > v(mm{£ y ,£ f )) - v(i y ) . (3) 

Observe that as X C F, it holds that £ x > £ y , so v{£ y ) > v{£ x ) as v is monotone decreasing. 
Now, one of the following holds: 

• ti£f<£y< 4 then Eq. © becomes u(£ f ) - i/(4) > u(£ f ) - u(£ y ), which holds. 

• If £ y < £f < 4 then Eq. (j3J) becomes v(£f) — v{£ x ) > ^(4) — v {^v)i which is equivalent 
to v{£f) > i/(4)- This in turn holds by the decreasing monotonicity of v. 

• If £ y < 4 < if then Eq. © becomes = z/(4) - u(£ x ) > u(£ y ) - v{£ y ) = 0. 

We conclude that Qp(-) is submodular. ■ 

To solve our problem, using the framework of Chekuri et al. |CVZ11] . we need the fol- 
lowing: 

(A) The target function is indeed submodular and can be computed efficiently This is 
Lemma 14.81 
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(B) State an LP that solves the fractional problem (and its polytope contains the optimal 
integral solution). This is just the original LP, see Eq. flTJ. 

(C) Observe that our rounding (i.e., contention resolution) scheme is still applicable in this 
case. This follows by Lemma 14.31 

One can now plug this into the algorithm of Chekuri et al. |CVZ11] and get an Q(a)- 
approximation algorithm, where a is the rounding scheme gap. The algorithm of Chekuri 
et al. |CVZllj uses a continuous optimization to find the maximum of a multi-linear exten- 
sion of the target function inside the feasible polytope, and then uses this fractional value 
with the rounding scheme to get the desired approximation. 
We thus get the following. 

Problem 4.9 Let P be a set of n points in the plane, and let F be a set of m objects in the 
plane. Let u(r) be a monotone decreasing function, which returns the amount of coverage a 
point gets if it is at distance r from one of the regions of F. Consider the scoring function 
that for an independent set H C F returns the total coverage it provides; that is, 

op(H) = 5>(<*p(H)) ■ 

P eP 

We refer to the problem of computing the independent set maximizing this function as the 
partial coverage problem. 

Theorem 4.10 Given a set of n points in the plane, and set m of unweighted objects in the 
plane with union complexity 0(u(n)), one can compute, in polynomial time, an independent 
set. Furthermore, this independent set provides an Q(n/u(n)) -approximation to the optimal 
solution of the partial coverage problem. 

Observe that the above algorithm applies for any pricing function that is submodular. In 
particular, one can easily encode into this function weights for the ranges, or other similar 
considerations. 

5 Weighted Rectangles 
5.1 The algorithm 

For the (original) independent set problem in the case of weighted axis-aligned rectangles, 
we can solve the same LP, where the set V (F) contains both intersection points and corners 
of the given rectangles. 

Define two subgraphs Si and S2 of the intersection graph: if the boundaries of fi and fj 
intersect zero or two times, put fcfj in Si; if the boundaries intersect four times instead, put 
fifj in 82- 

We first extract an independent set I of Si using the algorithm of Theorem 14.41 
It is well known (e.g., see |AG60j ) that S2 forms a perfect graph (specifically, a com- 
parability graph), so find a A-coloring of the rectangles of / in S2, where A denotes the 



1(3 



maximum clique size, i.e., the maximum depth in ,A(R). Let I' be the color subclass of / of 
the largest total weight. Clearly, the objects in I' are independent, and we output this set. 



5.2 Analysis 



As in Lemma 14 .1\ let H be any subset of F. Observe that if fifj £ 9i, then contains a 
corner of fj or vice versa. Letting Vj denote the corners of fj, we have 

x * x i ^EE XiX i ^mi x J x ^HH x i = Y. Ax ^ 

fi,fjeH,fifjeSi fjeHpeVj fieH,pefi /,? H P < r, /,? n.|> f /.•• H|>- \ : , / i i 

= 4£(H). 

Applying the same analysis as before, we conclude that the expected total weight of / 
(i.e., the independent set of Si) computed by the algorithm is of size fi(Opt). 
To analyze I', we need a new lemma which bounds the maximum depth of R: 

Lemma 5.1 A = O (log nj log log n) with probability at least 1 — 1/n. 

Proof: Fix a parameter t > 1. Fix a point p £ V (F). The depth of p in .A(R), denoted by 
depth (p, R), is a sum of independent 0-1 random variables with overall mean \i = 
1. By the Chernoff bound |MR95l page 68], 

Pr[depth(p, R) > (1 + 5)fi] < 



Xi < 



1+6 



for any 5 > (possibly large). By setting 5 so that t = (1 + <5)/i, this probability becomes 
less than (e/t) 1 . Since |V(F) | = 0(n 2 ), the probability that A > t is at most 0((e/t)*rj- 2 ), 
which is at most 1/n by setting the value of i to be 0(logn/ log log n). m 

By construction of we know that 

m > xyi w ' - ^ e^ ^ ~ ° pt ■ 1a> * 



/i6/' 



/ie-r 



where 1^ denotes the indicator variable for event A. Witht = (log n/ log log n), we conclude 
that 



E 



> f2(log log n/ log n) E u»j 

> ^(loglogra/logn) ■ Opt. 



Opt/n 
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5.3 Remarks 



Derandomization. This algorithm can also be derandomized by the method of conditional 
expectations. The trick is to consider the following random variable 



Opt • £ (! 

peV(F) 



r \depth(p,R) — t 
°p) ) 



where 5 P is the S from the proof of Lemma 15.11 and t is the same as before. This variable 
Z lower-bounds Ylf w i (^ ne bound is trivially true if A > t, since Z would be negative). 
Our analysis still implies that ~E[Z] > Q(log logn/ logn) ■ Opt (since the standard proof of 
the Chernoff bound |MR95l pages 68-69] actually shows that 



Pr 



depth(p, R) > (1 + S)fM 



E [(l + 5)depth(p,R)j 

< (1 + 5)^+^ - 



1 f 



S) 



1+5 



which for 5 = S p and fi < 1 implies E[(l + 5 p ) depth ( p ' R ) -< ] < (e/£)*). The advantage of work- 
ing with Z is that we can calculate E[Z] exactly in polynomial time, even when conditioned 
to the events that some objects are known to be in or not in R (since depth(p, R) is a sum of 
independent 0-1 random variables, making (1 + <5 p ) de P th (P' R ) a product of independent random 
variables) . 

We have thus proved: 



Theorem 5.2 Given a set of n weighted axis-aligned boxes in the plane, one can compute 
in polynomial time an independent set of total weight fi(log logn/ logn) ■ opt, where opt is 
the maximum weight over all independent sets of the objects. 



Higher dimensions. By a standard divide-and-conquer method |AvKS98] , we get an 
approximation factor of 0(log d_1 n/ log log n) for weighted axis-aligned boxes in any constant 
dimension d. 
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